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ABSTRACT 
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' Context, Far-infrared excess emission from main-sequence stars is due to dust produced by orbiting minor bodies. In these disks, larger 

CZ2 , bodies, such as planets, may also be present and the understanding of their incidence and influence currently presents a challenge. 

Aims. Only very few solar-type stars exhibiting an infrared excess and harbouring planets are known to date. Indeed, merely a single 
case of a star-planet-disk system has previously been detected at submillimeter (submm) wavelengths. Consequently, one of our aims 
is to understand the reasons for these poor statistics, i.e., whether these results reflected the composition and/or the physics of the 
planetary disks or were simply due to observational bias and selection effects. Finding more examples would be very significant. 
Methods. The selected target, q 1 Eri, is a solar-type star, which was known to possess a planet, q 1 Eri b, and to exhibit excess emission 
at IRAS wavelengths, but had remained undetected in the millimeter regime. Therefore, submm flux densities would be needed to 
better constrain the physical characteristics of the planetary disk. Consequently, we performed submm imaging observations of q 1 Eri. 
Results. The detected dust toward q 1 Eri at 870//m exhibits the remarkable fact that the entire SED, from the IR to mm-wavelengths, 
is fit by a single-temperature blackbody function (60 K). This would imply that the emitting regions are confined to a narrow region 
Cf") ' (ring) at radial distances much larger than the orbital distance of q 1 Eri b, and that the emitting particles are considerably larger than 

some hundred micron. However, the 870 yum source is extended, with a full-width-half-maximum of roughly 600 AU. Therefore, a 
physically more compelling model also invokes a belt of cold dust (17 K), located at 300 AU from the star and about 60 AU wide. 
Conclusions. The minimum mass of 0.04 M e (3 Mmooii) of 1 mm-size icy ring-particles is considerable, given the stellar age of J 1 Gyr. 
These big grains form an inner edge at about 25 AU, which may suggest the presence of an unseen outer planet (q 1 Eri c). 
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1 . Introduction belts and the zodiacal dust cloud in the solar system dMann et al.l 

[2006b . 

During the end stages of early stellar evolution, dusty debris 

disks are believed to be descendents of gas-rich protoplanetary For solar-type stars on the main-sequence, which are known 
disks. These had been successful to varying degrees in build- to exhibit infrared excess due to dust disks, one might ex- 
ing a planetary system. What exactly determines the upper cut- pect, therefore, a relatively high incidence of planetary sys- 
off mass of the bodies in individual systems, and on what time terns around them. Surveying nearly 50 FGK star s with known 
scales, is not precisely known. However, the presence of debris planet s for excess emission at 24 fjm and 70/iiri, iTrilling et all 

d2008h . 



around matured stars is testimony to the action of orbiting bod- d2008l) detected about 10-20% at 70jum, but essentially none at 

ies, where a large number of smaller ones are producing the dust 24jum, implying that these planetary disks are cool (< 100 K) 

through collisional processes and where a small number of big- and large (> 10 AU). However, in general, the conjecture that 

ger bodies, if any, are determining the topology (disks, rings and the infrared excess arises from disks lacks as yet observational 

belts, clumps) through gravitational interaction. The time evo- confirmation due to insufficient spatial resolution. In fact, until 

lution of the finer debris is believed to be largely controlled by very recently, there was only one main-sequence system known 

non-gravitational forces, though. By analogy, many debris disks that has an extended, resolved disk/belt structure and (at least) 

are qualitatively not very different from the asteroid and Kuiper one giant planet, viz. eEri, a solar-type sta r at the distance of 

only three parsec dGreaves et al.l ll998, 2005). Its planetary com- 

panion, e Eri b, has been detected indirectly by ast rometric and 

Send offprint requests to: R. Liseau radial velocity (RV) metho ds applied to the star (lHatzes et al.l 

* Based on observations with APEX, Llano Chajnantor, Chile. 120001 : iBenedict et al.l 120061) . whereas attempts to directly de- 
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tect the planet hav e so far been unsuccessful (lltoh et al.ll2006t 
Uanson et al.ll2007l) . 

As its name indicates, the object of the present study, q 1 Eri, 
happens to belong to the same celestial constellation of Eridanus, 
albeit at a larger distance (D = 17.35 ± 0.2 pc) and is, as such, 
unrelated to e Eri. The planet was discovered with t he RV tech- 
nique (for a recent overview, see iButler et"ai1l2006l) . These RV 
data suggest that the semimajor axis of the Jupiter-mass planet 
q 1 Erib is about 2 AU (Table[T|i. It seems likely that regions in- 
side this orbital distance have been largely cleared by the planet, 
whereas outside the planetary orbit, substantial amounts of ma- 
terial might still be present. 

In fact, IRAS and ISO data were suggestive of significant 
excess radiation above th e photospheric emission at w avelengths 
longward of about 20 /Ym. lZuckerman & Sondd2004l) interpreted 
these data in terms of dust in a disk at the orbital distance of 
30 AU and at a temperature of about 55 K. IChen et al.1 d2006l) 
fitted the far-infrared emission w ith the correspondin g values 
of 20 AU and 70 K, respectively. iTrilling et all (120081) derived 
20 AU and 60 K. In their entire sample of more than 200 stars, 
q 1 Eri (=HD 10647) has by far the highest 70 ^m excess. 

At mm-wavelengfhs. lSchutz et~aU d2005l) failed to detect the 
disk and assigned an upper limit to the dust mass of 6Mmooiv 
This is unsatisfactory, as the proper characterization of the dust 
around q 1 Eri would require valid long wavelength data. In the 
following, observations of q 1 Eri at 870 fim are described and 
their implications discussed. 

2. Observations and Data Reductions 

APEX, the Atacama Pathfinder Experiment, is a 12 m diameter 
submillimeter telescope situated at an altitude of 5100m on the 
Llano Chajnantor in northern Chile. The telescope is operated 
by the Onsala Space Observatory, the Max-Planck-Institut fur 
Radioastronomie, and the European Southern Observatory. 

The La r ge A pex BOlometer CAmera (LABOCA, 
ISiringo et al.l [20071) is a multi-channel bolometer array for 
continuum observations with 60 GHz band width and centered 
on the wavelength of 870yt/m. The array, having a total field 
of view of 11', is spatially undersampled and we therefore 
adopted spira l patte rn observing as the appropriate technique 
dSiringo et al.l 120071) . This procedure results in fully-sampled 
maps with a uniform noise distribution over an area of about 
8'. During the nights of August 1-4, 2007, we obtained 32 such 
individual maps, for about 7.5 min each with central coordi- 
nates RA=01 h 42 m 29 s 32 and Dec=-53° 44' 27'.'0 (J2000). The 
LABOCA beam width at half power (HPBW) is 18"6 ± l'.'O. We 
focussed LABOCA on the planet Jupiter and the rms-pointing 
accuracy of the telescope was 3" to 4". 

W e reduced the data with the BoA software dSiringo et alj 
|2007|) . which included flat fielding, baseline removal, despiking 
and iteratively removing the sky noise, and filtering out the low 
frequencies of the 1 //-noise, with the cut-off frequency corre- 
sponding to several arcminutes. The software also accounts for 
the map reconstruction and the absolute calibration, using the 
opacities determined from numerous skydips (zenith opacities 
were in the range 0.1 to 0.3) and observations of the planets 
Uranus and Mars. The final result is an rms-noise-weighted av- 
erage map (Fig.Q]). 



3. Results 

The final product of the reduction process is the 870//m im- 
age presented in Fig.Q] which shows the central 5' x 5' of 
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Fig. 1. q 1 Eri observed at 870//m with the submm camera 
LABOCA at the APEX telescope (HPBW- 18"). Within the 
positional accuracy, the star is at the origin of the image (see 
TableEl referring to 01 h 42 m 29 s 32, -53° 44' 27"0, J2000.0) and 
offsets are in arcsec. The lowest contour corresponds to 2cr of the 
rms noise of the flux density, F v , and increments are in steps of 
lcr. The colour coding, in units of Jy/beam, is shown by the scale 
bar to the right of the image, which has been smoothed with a 
circular 18" Gaussian. At the distance of the star, the space be- 
tween two tick marks (=20") corresponds to 350 AU. 

Table 1. Physical properties of the star and its planet * 



Parameter 



Value 



The star q 1 Eri 

Distance, D 

Spectral type and luminosity class 

Effective temperature, T eff 

Luminosity, L 8ta r 

Surface gravity, log g 

Radius, /? star 

Mass M star 

Metallicity, [Fe/H] 

Age 

The planet q 1 Erib 

Period, P 

Semimajor axis, a or bii 
Eccentricity, e 
Mass, M sin i 



17.35 pc 
F8-9V 
6100K 
1.2 Lo 

4.4 (in cm s~ 2 ) 
1.1 

1.1 Mo 
-0.08 

(> l-2)Gyr 



2.75 ± 0.15 yr 
2.0±0.2AU 
0.2 ± 0.2 
0.9 + 0.2M Jupta . 



See lButler et al.l d2006h and references cited in the text. 



the LABOCA map. The peak flux in the map is found at the 
position of q 1 Eri and a few other pointlike features of low 
intensity are also present, one of which is close to the star. 
If not merely noise, these low signals could be due to extra- 
galactic background sources, as the displayed n umber density 
is consistent with that observed elsewhere (e. g., lLagache et al.l 
120051: iBertoldi etaT] 120071: llvison et alj|2007l) . Other, comple- 
mentary observations (e.g., optical, IR, X-rays, radio interfer- 
ometry) would be required for their identification. 
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Table 2. Physical properties of the q 1 Eri dust system 



Parameter 


Value 


r CoR UllScL , [LlLi, L\0 ) 


, tj j, terror. ) 


P^qV fliiv H^ncif\r n F (C\ CY\ ) — 870 //m 

rcdK. iiua ueiisiLy , r v \\j, vt), a — o/UjUiii 


^ i o.z it \j.o / [iu yi ueaiii 


Integrated flux density, j F v dadS 


±4.i^ mjy, t v > z<x 


Major axis" (FWHM) 


1,1" ±2" (640 + 35)AU 


Position angle", pa 


55° ± 4° (north over east) 


Minor axis" (FWHM) 


23" ± 1" 


Inclination angle, i 


> 52° (90° = edge-on) 


Fractional luminosity, L bb /L star 


1.1 x 10- 4 


Inner (outer)* Temperature, 


60K(17K) 


Inner (outer) 6 Radius, r b b 


25 AU(300AU) 


Inner (outer) 6 Width, Ar bb 


0.02AU(60AU) 


Inner (outer) 6 Minimum mass'', M dllst 


0.04 M e (0.15 M ) 



" Two-dimensional Gaussian fits with la formal fitting uncertainties. 
6 An outer dust belt is implied by the extent of q 1 Eri at 870 pm. 

2cm 2 g~' (p = 1.18 gcirT 3 , a milx = 1 mm, n{a) cc a -3 - 5 ). 



fl0 115 Hz 



The derived flux densities of q 1 Eri are provided in Tablef2j 
which presents the results from fitting the data to a two- 
dimensional Gaussian function. The indicated errors are formal 
fit errors only, based on lcr rms values. In individual cases, e.g., 
pa, realistic errors could be twice as large. The error on the in- 
tegrated flux density in Tableff] also includes an uncertainty of 
10% in the absolute calibration. 

The 870 fim source is at best only marginally resolved in the 
North-South direction (formal fit result is 23" + 1"), whereas 
it is clearly elongated in approximately the East- West direction 
(37" + 2"). At the distance of q 1 Eri, this corresponds to a disk 
diameter of 640 AU and assuming a circular shape, these disk 
dimensions yield an inclination with respect to the line of sight 
of i > 52°, not excluding the possibility that the disk is seen 
essentially edge-on (i close to 90°). The vertical disk scale height 
is undetermined. 



4. Discussion 

4.1. Physical conditions and the age of the system 
The spectral type of q 1 Eri is slightly earlier than that of 



Ihe spectral type ot q is slightly e arlier than that or 
the S un (F8-9 V. iNordstrom et all l2004t iDecin et"aT] l2000l 



120031 IZuckerman & Song|l2004l Che n et al. l l2Q06hT with the ef- 
fective temperature bei ng bracketed by the extremes 6040 K 
(INordstrom et alj|2004l) and 6260 K dChen et al.l 120061) , with 
the mean of 6150 K, i.e., essententially the value given by 
iButler et al.ld2006l 6105 K) (see Table©. 

Literature estimations of likely ages for q 1 Eri span the 
range 0.3 to 4.8 Gyr (with an entire range of 0.0 to 7.0 Gyr , 
Decin et al.ll2000t IZuckerman & Sonill2004t IDecin et alj|2003t 



Nordstrom et alj|2004t IChen et al.ll2006l) . However, the level of 



chromospheric activity (logR' - -4.7) suggests an age of 
1.9 Gyr (see Eq. 15 of I Wright etaTl |2004l) . The star has also 
been detected in X-rays with ROSAT (logL x = 28.3, J . Sanz , 
private communicatio n), yielding 1 .2 Gyr dRibas et al.l 120051: 
iGuinan & Engiel 120071) . This value is also consistent with the 
stellar rotation period of about 10 days (uncorrected for sin/; 
lEcuvillon et al.l 120071) . It is clear that the star is definitely on 
the main-sequence and that the age of the system likely exceeds 
10 9 yr. 
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Fig. 2. The fit to the SED of q 1 Eri (HD 10647) is shown from 
1000 A to beyond 1 mm in red. The stellar photosphere is repre- 
sented by a Kurucz ATLAS 9 model atmosphere and the excess 
emission by a single temperature blackbody curve (in black). 
Also, shown in red, are Simbad, TI MMI2 11.9 ^m, a nd IRAS 
FSC data, in addition to ISO 60 fim dDecin et alj|2000h . Spitzer 
data have been given higher w eight and are sho wn in green, 
viz. 8.5-13 yum and 30-34 jum dChen et al.l 120061) . and 24yum 
and 70 fim dTrilling et al.ll2008l). respectively. T he SEST-SIMBA 
3cr upper limit at 1 .2 mrn ifSchiitz et"aT1l2005h and our APEX- 
LABOCA point at 870 /an are shown in blue. The black dots 
refer to the combined 60 K modified blackbody (J3 = 1) and 
17 K blackbody (see the text). 



4.2. The nature of the emitting particles 

The absence of spectral features in the 10 to 30 fim region sug- 
gests that the dust grains are consid erably larger than 10//m 
dChen et alj 120061: ISchiitz et all 120051) . Remarkably, the spec- 
tral energy distribution (SED) of the excess emission can be 
fit by a single-temperature blackbody of 60 K, from the in- 
frared to the submm/mm regime (see Fig.|2|. The blackbody 
character is determined by the LABOCA flux and independent 
of the relative weights assigned to the mid- and far-infrared 
data. The radial distance from the central star, at which a grain 
has attained thermal equilibrium, is approximately given by 
[(1 -A)/(16?r e&) (L star /T 4 )] 1/2 , where A and e are the integrated 
reflectivity and emissivity, respectively. For a blackbody this re- 
duces to rbb = (^star/2) (r e ff/Tbb) 2 , which for 7bb = 60 K, yields 
a minimum distance of 25 AU for the q 1 Eri dust (Table|2|. Taken 
at face value, this would mean that the range in dust tempera- 
tures is very limited: single values lower than 50 K or as high 
as 100 K can be excluded (without giving higher weight to the 
Spitzer data, Tbb becomes closer to 70 K). Therefore, rbb, is de- 
termined to better than within a factor of two. For unit filling 
factor, the blackbody emitting regions would appear to be con- 
fined to a very narrow ring-like structure (see Table|2]). 

The fractional luminosity is 1 x 10~ 4 and the emission is 
optically thin. The blackbody fit also implies that the emitting 
particles have sizes largely in excess of 100 /im (27ra > A) and 
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that these grains have grey opacities in the infrared to submm, 
i.e., k + k(A). Given the available evidence, it is not possible, 
however, to tell the actual sizes of the particles or their absolute 
opacities. 

Some insight m a y be gained from the work of 
Miyake & Nakagawa] d 19931) . who explored the optical 
properties of dust that produces small values of the opacity 
index, and presented opacities over a broad range in frequency 
and particle size. Maximum opacity, max/c ~ 2cm 2 g _1 , was 
found for the size a max = 1 mm at v = 10 1L5 Hz (A ~ 1 mm) 
and for larger particles, k decreases rapidly (as ~ 1/ ya max ). 
This assumes compact spheres of density p = 1.18 gem -3 
(well-mixed silicates and water ice) and being distributed in 
size according to n{a) oc a 1 ', with p = -3.5. The adopted 
density is consistent with values determined for Kuiper Belt 
objects dGrundv et al.ll2007l and references therein). The value 
of K\ mm is not strongly depen dent on p, as long as -4 < p < -2 
dMivake & Nakagawal 119931) . I n general, these results are in 
agreement with other w ork (e.g jKrugel & Siebenmorgenll994t 
IStognienkoetalJJ.19951) . 

For this maximum value of k, the 870jum flux density 
yields a minimum mass Md us t = F v D 2 /K v B v (Td ust ) > 3 M M „ n 
(0.04 M e , see Table[2]i. This minimum mass is larger than the 
'blackbody mass' one would infer from the effective area of the 
blackbody (6.85 x 10 26 cm 2 ) and f or the same a and p, but con- 
sistent with the (re-scaled) result of lSchiitz et al.l d2005l) . 

The fact that the 870/im source appears linearly resolved, 
speaks against the narrow ring scenario, and the existence of a 
more massive and colder belt (T — 17 K at r — 300 AU, say) 
can at present not be excluded (see Fig.|2|. The relative width of 
such a cold belt would seem less implausible, viz. Ar/r > 0.2 
and, hence, its physical width would be at least 60 AU. Also, 
the equilibrium temperature of the dust would be esssentially 
constant. With the same parameters as before, the mass would 
scale simply as the ratio of the temperatures, yielding 13 Mmooiv 
Of course, at shorter wavelengths, the spectrum would have to 
be steeper than the blackbody SED, implying values of /3 > 0, 
where /3 parameterizes the frequency dependence of the opacity, 
i.e., k v oc v 13 . Together with the 17 K blackbody, the SED can 
be fit with a modified 60 K blackbody with = 1 .0 longward of 
100 pm (Fig. [2]). This model would physically be more attractive, 
but by its ad hoc nature would be of course not unique, and bet- 
ter constraints on the physical parameters would require a better 
sampling of the SED. 

4.3. Relation to the planet q 1 Eri b 

An in dex of the order of -3.5, used by iMiyake & N akagawa 
d 1993b . could indicate that the size distri bution resulted from a 
collis ional cascade (for a discussion, see iThebault & Augereaul 
120071) . The observed absence of small debris (a of the order of 
1 pm or smaller) in the q 1 Eri disk suggests that its production 
has ceased and that it had diminished on a short time scale com- 
pared t o the age of th e system. Remarkably, the density param- 
eter of IWvattl d2005l) has a value that, given the age of q 1 Eri, 
is atypically large (770 ~ 1000 for the correct form of Eq. 7), 
which may mean that the observed absence of warm material 
is caused by radiation pressure blow-out, rather than by the ac- 
tion of a planet, hindering the migration inward toward the star. 
Anyway, at 2 AU distance, the Jupiter-size planet q 1 Erib could 
hardly have had any influence on particles o ut to the inne rmost 
edge at 25 AU and on orbits far beyond that dWvattll2005l) . The 
existence of this belt of large grains may point to the presence 



of another major planet. It would therefore seem important to 
verify or to disprove the existence of q 1 Eri c. 

5. Conclusions 

Below, our main conclusions from this work are summarized: 

• Observations of the solar-type star q 1 Eri and its planet 
q 1 Erib at 870//m revealed a source with peak emission at 
the position of the star. The source appears extended to the 
LABOCA beam, i.e., elongated in roughly the East- West di- 
rection (640 AU) but essentially unresolved in the perpendic- 
ular direction. 

• At an age exceeding 1 Gyr, the fractional luminosity of the 
infrared excess is very high (> 10~ 4 ). The entire SED of this 
excess emission, extending from J 20 pm to 1 mm, is fit by a 
single-temperature blackbody (7bb = 60 K). 

• This would imply that the emitting regions are located at 
about 25 AU from the star and in addition, very limited 
in spatial extent (ring-like). Exhibiting a grey opacity over 
the entire wavelength range, the emitting particles must be 
large (>> 100 pm). Using the theoretically derived maxi- 
mum grain opacity for 1 mm-size icy particles, we estimate 
a minimum mass of the dust belt of 0.04 M m . 

• It seems highly unlikely that the planet q'Erib at 2AU 
would be responsible for the clearing of the region from 
small dust particles interior to 25 AU, but may hint at the 
existence of another planet. 

• Taking the observed extent of the 870 pm source into account 
leads to an emission model in which an outer cold (17 K) 
dust belt needs to be included. This belt would be centered 
on the radial distance of 300 AU and have a width of 60 AU. 
This belt is inclined at ; > 52°, possibly viewed at an angle 
close to edge-on. 

• The extreme character of the debris disk around the relatively 
old star-planet system q 1 Eri provides yet another example of 
the large diversity of such disks. 

Acknowledgements. We are indebted to the staff at the APEX facility in Sequitor 
and at the telescope site for their enthusiastic and skillful support during our 
observing run. The thoughtful comments by the referee were much appreciated. 
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